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£h ■ Abstract 

We study the classical and quantum models of a Friedmann-Robertson- Walker (FRW) cos- 
mology in the framework of the gravity theory proposed by Hofava, the so-called Hofava-Lifshitz 
theory of gravity. Beginning with the ADM representation of the action corresponding to this 
qj', model, we construct the Lagrangian in terms of the minisuperspace variables and show that in 

comparison with the usual Einstein-Hilbcrt gravity, there are some correction terms coming from 
^' ( , the Hofava theory. Either in the matter free or in the case when the considered universe is filled 

bi)' with a perfect fluid, the exact solutions to the classical field equations are obtained for the flat, 

closed and open FRW model and some discussions about their possible singularities are presented. 
We then deal with the quantization of the model in the context of the Wheeler-DeWitt approach 
of quantum cosmology to find the cosmological wave function. We use the resulting wave functions 
0^ ' to investigate the possibility of the avoidance of classical singularities due to quantum effects. 
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1 Introduction 

Since 2009 when a new theory of gravity was first introduced by Hofava [1], many efforts have 
been made in this area and the corresponding results have been followed by a number of works, the 
main motivations of which lie in the results of quantum gravity and cosmology, see for instance [2] . 
Since this kind of gravitational theory has its roots in the Lifshitz work on the second-order phase 
transition in solid state physics, it is usually referred to as the Hofava-Lifshitz (HL) theory of gravity. 
Like another candidates for quantum gravity such as string theory, HL gravity is also a completion 
of general relativity (GR) at high energy ultraviolet (UV) regime and reduces to standard GR in 
the low energy infra-red (IR) limit. However, in the framework of the HL theory, the well-known 
phenomenon of Lorentz symmetry breaking at high energies is described somehow in a different way. 
Indeed the basic idea behind HL is that the Lorentz symmetry will be broken through a Lifshitz-like 
process, i.e., through an anisotropic scaling (characterized by a scale parameter b and the dynamical 
critical exponent z) between space and time as 

t -> b z t, x -4 bx. (1) 

There are theories correspond to the different values of z. While for z = 1 the standard relativistic 
scale invariance with Lorentz symmetry (the IR limit) is recovered, the UV gravitational theory 
proposed in [1] requires z = 3. Because of the asymmetry of space and time in HL theory the most 
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common representation of the space-time metric is its ADM form which in terms of the lapse function 
iV(t,x), shift vector iV a (i, x) and spacial metric h ab (t,x.) takes the general form 

-iV 2 (t,x) + iV a (t,x)iV a (t,x) N b (t,x) \ ( , 

JV a (t,x) h ab (t, X ) J' { > 

Depending on whether the lapse function is a function only of t or of (t,x), the theory is called 
projectable or non-projectable. As pointed out in [1], since in cosmological metrics the lapse function 
labels the time parameter, it is quite reasonable to choose it as a projectable function. However, 
more general cases in which the lapse function is taken as a non-projectable function are studied 
in [3]. There is yet another issue related to the HL theory, the so called detailed balance condition. 
The gravitational action of the model consists of the kinetic part Sk and the potential part Sy as 
S = Sk + Sy- The kinetic part comes from the Einstein-Hilbert action usually written in terms of 
the ADM variables (see the relation (11) bellow). For the potential part the following general form 
is proposed in [1] 

Sv = f d 4 x^V[h ab ], (3) 

where V is a scalar function which depends only on the spacial metric h ab and its spacial derivatives. 
Taking a three dimensional HL model with z = 3, this function may be constructed by a superposition 
of the quadratic (in curvature) terms such as 

V a R bc V a R bc , V a R bc V b R ca , RV 2 R, R ab V 2 R ab , (4) 

and cubic terms such as 

p3 pa p6 pc pp nab /r\ 

K , K b K c K a , KK ab K . (b) 

Among the very different possible combinations, Hofava considered a special form, known as "detailed 
balance condition" satisfying model as 

S v ~ J d 4 x^E ab G abcd E cd , (6) 

where Gabcd is the DeWitt metric and E ab is a tensor constructed by variation of some function W[/i &] 
with respect to the spacial metric. By choosing a suitable ansatz for the function W, he then showed 
that in a z = 3 theory with detailed balance condition the potential is a combination of the terms 
(see [1] for details) 

V a R bc V a R bc , V a R bc V b R ac , V a RV a R. (7) 

The issue of the detailed balance condition makes the theory to have simpler renormalization prop- 
erties. Although the detailed balanced system exhibit simpler quantum behavior, it is shown in [4] 
that if one relaxes it, the resulting Lagrangian with extra allowed terms is well-behavior enough to 
recover the model with detailed balance. 

In this paper we shall consider a FRW cosmological model in the framework of a projectable HL 
gravity without detailed balance condition. Our approach to deal with such a problem is through 
its representation with minisuperspace variables. Minisuperspace formulation of HL cosmology is 
studied in some works, see for instance [5]- [7], to obtain its possible classical and quantum solutions. 
Here, we first consider the vacuum case in which the cosmological model is free from the presence 
of any kind of matter fields and see that the corresponding classical solutions exhibit some types of 
singularities. We then will add a perfect fluid as the matter into the model and obtain some exact 
solutions in the cases of flat, closed and open FRW cosmologies. Since our aim in the quantum part 
of the model is to investigate the time evolution of the wave function, we prefer to use the perfect 
fluid in its Schutz formalism in which the Hamiltonian of the fluid consists of a linear momentum, 
the variable canonically conjugate to which may play the role of a time parameter. In both vacuum 
and perfect fluid cases, we construct the corresponding quantum cosmology based on the canonical 
approach of Wheeler-DeWitt theory to see how things may change their behavior if the quantum 
mechanical considerations come into the model. 
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2 The model 



In this section we start by the FRW cosmology within the framework of HL gravity. In a quasi- 
spherical polar coordinate the metric of space time is assumed to be 



ds z 



-N 2 (t)dt 2 + a 2 (t) 



dr 2 



1 — kr 2 



(d'd 2 + sin 2 &d<p) 



(8) 



where N(t) is the lapse function, a(t) the scale factor and k = 1, and —1 corresponds to the closed, 
flat and open universe respectively. It is clear that in terms of the ADM variables the above metric 
can be written as 

ds 2 = -N 2 (t)dt 2 + h ab dx a dx b , (9) 



in which 



h ah = a 2 (t)diag 



1 



; ,r 2 ,r 2 sin 2 i? I , 



1 — kr 2 

is the intrinsic metric induced on the spatial 3-dimensional hyper surf aces. The action of the model 
consists of the gravitational part S g and the matter action S m as 



(10) 



The matter part of the action is independent of the HL corrections to the gravity part. In the context 
of the ADM formalism, following [4] and with the same notation as is used in [5, 6], the projectable 
HL gravity without detailed balance has the (gravitational) action (in what follows we work in units 
where c = h = I) 



M 2 Pl 



[ d 3 xdtNVh \K ab K ab - XK 2 - g M 2 Pl - 9l R - M p 2 (g 2 R 2 + g 3 R ab R ab 
JM 1 v 



-Mpf (g A R 3 + g 5 RR a b R b a + g 6 R a b R b c R c a + g 7 RV 2 R + g s V a R bc V a R bc ) 



+ M 2 Pl [ d 3 x^hK, (11) 



where Mpi = (87rG) -1 / 2 is the Planck mass, K ab are the components of the extrinsic curvature tensor 
which describes how much the spatial space h ab (which is the boundary dM of the four-dimensional 
manifold J\A) is curved in the way it sits in the space time manifold. Also, h and R are the determinant 
and Ricci scalar of the spatial geometry h ab respectively, and K represents the trace of K ab . The 
constants A and gi (i = 0, 1, ...,8) denote the HL corrections to the usual GR. In comparison with 
the ADM representation of GR action, that is 



Mpi 



f d 3 xdtNVh (K ab K ab - K 2 + R - 2A) + 2 / d 3 xVhK 
JM V ' JdM 



we conclude that the consensus GR is recovered if one sets the cosmological constant as A = goM Pl /2, 
<7i = — 1 and A = 1. However, as is mentioned in [5], we take A as a running constant which represents 
the IR limit of the gravitational theory. With these identifications, the action (11) may be rewritten 

as 



c _ M Pl 



J d 3 xdtN^h [K ab K ab - XK 2 + R-2A- M p 2 (g 2 R 2 + g 3 R ab R ab ^j 



-M p ? [g 4 R 3 + g 5 RR a b R b a + g Q R a b R b c R c a + g 7 RV 2 R + g s V a R bc V a R bc ) 
From its standard definition, the extrinsic curvature is given by 



+ M 2 Pl 



dM 



d 3 xVhK. (12) 



K, 



ab 



1 

2N 



N alb + N bla -^) 



(13) 
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where N a is the shift vector and N a \ b represents the covariant derivative with respect to h^. Since the 
shift vector is absent in the FRW models, a simple calculation based on the above definition results 
in K a i)K ab = -jj^jz and K = — jj^, where a dot represents differentiation with respect to t. Also, the 
Ricci tensor and the Ricci scalar correspond to the 3-geometry h a i, can be obtained as R a b 



2kh„ 



and R = 0%. The gravitational part for FRW model may now be written by substituting the above 



results into action (12), giving 



3V -Mp,(3A-l) 



J" tN \-jp 



+ 



6ka 



2Aa 3 



3(3A - 1) 3(3A - 1) 



P l 



I2k 2 (3g 2 + g 3 ) 



-M 



-4 
PI 



3a(3A - 1) 
24k(9g 4 + 3g 5 + g & \ 



where Vq = J d 3 xVh is the integral over spatial dimensions. If we set 
led to the point-like Lagrangian 



3a 3 (3A - 1) 

■jV M 2 pl (A\-l) 
2 



(14) 



C g = N 



aa 



+ g c ka-g A a - 
N z a 



g r k 2 g s k 



1, then we are 
(15) 



where its coefficients are defined as [8] 
2 2A 



9c 



3A - 1 



9A 



3(3A - 1) 



g r = <5V (3g 2 + g 3 ), g s = 18V 2 (3X - 1)(9<? 4 + 3g 5 + g 6 ). (16) 



Now, the gravitational part of the Hamiltonian for this model can be obtained from its standard 
definition H g = aP a — C g . Noting that 



d£g 

da 



aa 
N' 



one gets 



H g = NT-Lg = N 



Pa i i 3 , 9rk 2 g s k 

— g c ka + g A a 3 H h -^r 

4a a a 6 



(17) 
(18) 



We see that the lapse function enters in the Hamiltonian as a Lagrange multiplier, as expected. Thus, 
when we vary the Hamiltonian with respect to N, we get H g = 0, which is called the Hamiltonian 
constraint. On the classical level this constraint is equivalent to the Friedmann equation, while 
on the quantum level, the operator version of this constraint annihilates the wave function of the 
corresponding universe, leading to the so-called Wheeler-DeWitt equation. Now, let us deal with the 
matter field with which the action of the model is augmented. As we have mentioned, the matter part 
of the action is independent of modifications due to the HL terms. Therefore, the matter may come 
into play in a common way and the total Hamiltonian can be made by adding the matter Hamiltonian 
to the gravitational part (18). To do this, we consider a perfect fluid whose pressure p is linked to its 
energy density p by the equation of state 

p = up, (19) 

where — 1 < uj < 1 is the equation of state parameter. According to Schutz's representation for the 
perfect fluid [9], its Hamiltonian can be viewed as (see [10] for details) 

Pt 



N- 



,3u) : 



(20) 



where T is a dynamical variable related to the thermodynamical parameters of the perfect fluid and 
Pt is its conjugate momentum. Finally, we are in a position in which can write the total Hamiltonian 
H = H„ + H m as 



H = NU a 



N 



p2 
± a 

' 4a 



j 3 Qr & 

g c ka + g A a + 



a 



1 <l± I'r 



a° 



a 



3u> 



(21) 
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Figure 1: Qualitative behavior of the scale factor for the flat universe. Solid and dashed lines denote the relation (23) 
for positive and negative signs respectively. 

The setup for constructing the phase space and writing the Lagrangian and Hamiltonian of the model 
is now complete. In the following sections, we shall deal with classical and quantum cosmologies which 
can be extracted from a theory with the above mentioned Hamiltonian. 

3 Cosmological dynamics without matter 
3.1 Classical model 

When the model is free of the contribution of any kind of matter fields, the dynamics is described 
by Hamiltonian (18). Therefore, the classical equations of motion is governed by the Hamiltonian 
equations. Equivalently, we can use the constraint equation % = which is nothing but the variation 
the action with respect to N . In this sense, from (17) and (18) we obtain 

ah 2 + g c ka - g A a 3 - - ^ = 0. (22) 
a a 6 

For the flat case k = 0, there is no difference between the HL cosmology and the usual FRW model 
except that the cosmological constant shifts as A — > 3g\ = 3 ^_i • In this case the above equation 
admits the following solutions 

a(t) = a e ± ^ t , (23) 

where a$ is an integration constant and the positive (negative) sign in the power of the exponential 
function corresponds to the expansion (contraction) universe. Thus, in this case we have two distinct 
branches of solutions, one of which begins with zero size at t = — oo and expands forever according to 
an exponential de Sitter law, while another has an opposite behavior, i.e., decreases its size from the 
large values at t = — oo and tends exponentially to zero as t — > oo. The situation of the dynamical 
behavior of the scale factor in this case is shown in figure 1. 

For k 0, equation (22) does not seem to have analytical solutions. So we consider the behavior 
of its solutions only for some limiting cases. When the scale factor is very small, we keep the two last 
terms in (22) and rewrite it as 

■;2 _ _ 

a a° 

from which we obtain the following implicit relation between t and a 



ah 2_9_r^_9_sk = ^ (24) 

"TV 



a^g r a 2 + kg s k 9s In (g r a + ^ g l r d z + kg. 



25r 2gT 



±t - t , (25) 



where to is an integration constant. In figure 2, we have plotted the scale factor versus time based 
on the above relation. 

On the other hand, for the late time of cosmic evolution which the scale factor is expected to be 
large the terms with coefficients g c and g\ in (22) become more important. In such a case we write 
this equation as 

a 2 + g c k - g A a 2 = 0. (26) 
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Figure 2: The figures show the evolutionary behavior of the non-flat universe in the early times based on (25). The 
left and right figures correspond to the k = 1 and k = — 1 respectively. 




Figure 3: Left: Late time behavior of the non-flat universe for k = 1. Right: The same figure for k = — 1. 



For k = +1, this equation has two sets of solutions as 



2c/a 



Sa(*+*o) 



and 



a n (t) 



1 



(27) 



(28) 



where to is an integration constant. Each of these solutions consists of two branches. In one branch the 
universe contracts and when reaches a minimum size undergoes to an expansion period. Therefore, 
in the case of of a closed universe we have bouncing cosmologies in which the bounce occurs at 



t -- 
k 



[to + 



„ , — ln( )1 for ar(t) and at t 

^■\fgK k 9c9a n 1 w 



\n(g c g\)] for ajj(t). On the other hand, for 



-1 the solutions to the equation (26) take the form 



a 1 it) 



a(t) 



l 

2</A 



Sa(*-*o) 



9A(t+t ) 



9c9Ae 



9c9\e 



9A(t-to) 



9A(t+to) 



t < t 



t>-t + 



ln(gc3A) 



(29) 



We see that in the time interval t < to- 



ln(g c gx) the scale factor decrease and tends its evolution 



at to 



■■ ln(g c gA) with zero size. Also, for t > —to + , ]n(g c gA), the universe has an expanding 



behavior begins its evolution with a zero size singularity at —to + 



\n(g c g\). These two branches 



of solutions are separated from each other by a classically forbidden region corresponding to the time 



interval to 



ln(<7 c <7A) < t < — to + 7= ln(g c g\) for which no open classical solutions exist. The 



above results are summarized in figure 3. 



3.2 Quantum model 

Now, we shall study the quantum behavior of the model described by the Hamiltonian (18). One way 
to do such a study is to investigate the well-known Wheeler-DeWitt equation for the corresponding 
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universe, that is 



a da 2 



p d 
a 2 da 



+ 4 -g c ka + g A a 3 + 



9rk^ g s k 



VP(a) = 0, 



(30) 



where ^(a) is the wave function corresponds to the quantum universe and the parameter p represents 
the ambiguity in the ordering of factors a and P a in the first term of (18). In the case of a flat universe 
k = 0, to see the correspondence of the classical and quantum solutions, we note that for large values 
of scale factor, the behavior of the system can be obtained in the Wentzel-Kramers-Brillouin (WKB) 
(semiclassical) approximation. Then substituting ^(a) = 0,(a)e lS ^ in equation (30) leads to the 
modified Hamilton-Jacobi equation 



1 

"4o 



Ta) +9Aa 3 + Q(a) 



H[a,P a 



0, 



(31) 



in which the quantum potential is defined as Q(a) = ■An'- 



ianlu? - Jn^- lt[s wel1 known that the 
quantum effects are important for small values of the scale factor and in the limit of the large scale 

factor can be neglected. Therefore, in the semiclassical approximation region we can omit the Q term 

.2 /— 3 



in equation (31) and obtain the phase function S(a) as S(a) = ±^^/gXa a . In the WKB method, 
the correlation between classical and quantum solutions is given by the relation P a = ^ . Thus, 
using the definition of P a in (17), the equation for the classical trajectories becomes a(t) = a^^*^ 1 , 
which shows that the classical cosmology of equation (23) is exactly recovered. The meaning of this 
result is that for large values of the scale factor, the effective action corresponding to the expanding 
and contracting universes is very large and the universe can be described classically. On the other 
hand, for small values of the scale factor we cannot neglect the quantum effects and the classical 
description breaks down. Since the WKB approximation is no longer valid in this regime, one should 
go beyond the semiclassical approximation. For the flat FRW metric we have k = and the two 
linearly independent solutions to equation (30) can be expressed in terms of the Bessel functions 
J v {x) leading to the following general solution 



*(a) 



c l J l/3 Q\/5A a3 ) +C2J-1/3 (J\/5A a " 



(32) 



where c\p are integration constants and we have taken p = 1. With an eye to the classical solutions 
(23), it is clear that the model is free of classical Big-Bang singularity in which a(to) = for some 
to- Therefore, we expect that the quantum model predicts a zero probability for the creation of 
the universe with zero size. This fact lead us to impose the boundary condition ^f(a = 0) = on 
the wave function, which results in C2 = 0. Note that equation (30) is a Schrodinger-like equation 
for a fictitious particle with zero energy moving in the field of the superpotential U(a) = —4gAa 4 . 
Usually, in the presence of such a potential, the mini-superspace can be divided into two regions, 
U > and U < 0, which could be termed the classically forbidden and classically allowed regions, 
respectively. In the classically forbidden region the behavior of the wave function is exponential, while 
in the classically allowed region the wave function behaves oscillatorily. In the quantum tunneling 
approach [11], the wave function is so constructed as to create a universe emerging from nothing by 
a tunneling procedure through a potential barrier in the sense of usual quantum mechanics. Now, 
in our model, the superpotential is always negative, which means that there is no possibility of 
tunneling anymore, since a zero energy system is always above the superpotential. In such a case, 
tunneling is no longer required, as classical evolution is possible. As a consequence, the wave function 
always exhibits oscillatory behavior. In figure 4, we have plotted the square of the wave functions 
for typical values of the parameters. It is seen from this figure that the wave function has a well- 
defined behavior near a = and describes a universe emerging out of nothing without any tunneling. 
Now to see that how the quantum solutions may describe an expanding or contracting universe, we 
use a mechanism which we have called the probabilistic evolutionary process (PEP), based on the 
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Figure 4: Qualitative behavior of the probability (= |\P(a)| 2 ) versus scale factor for the flat empty universe, see relation 
(32). 

probabilistic structure of quantum systems, to provide a sense of the evolution embedded in the wave 
function of the universe (see [12] for details). This is based on the fact that in quantum systems 
the square of a state defines the probability, V a = \^(a)\ 2 . To make the discussion more clear, let 
us take a specific initial condition corresponding to the point P. Then, PEP states that the system 
(here specified by the scale factor a) moves continuously to a state with higher probability and thus 
P moves to the right to reach the point Q, a local maximum. Therefore, the (expanding) universe 
begins its evolution by a monotonically increasing scale factor to reach the point Q. In this sense 
the transition R — >■ Q describes a contracting universe. On the other hand, the emergence of several 
peaks in the wave function may be interpreted as a representation of different quantum states that 
may communicate with each other through tunneling. This means that there are different possible 
universes (states) from which the present universe could have evolved and tunneled in the past, from 
one universe (state) to another. Based on this interpretation, one can argue that from the states 
located on the larger peaks of the wave function a tunneling may occurs to reach the states located 
on the smaller peaks and then the universe evolved according to PEP. We schematically show such 
procedures as 

,. . J-.PEP s^tunnel.r^PEPjj /nn\ 

expanding universe : P — >Q — > 1 — >U... (66) 

As the scale factor grows the probabilities become small and smaller and the universe undergoes its 
classical region. For the contraction universe similar discussion as above would be applicable as well. 

Now, let us to deal with the Wheeler-DeWitt equation (30) in the case of k ^ 0. Like the classical 
solutions in this case, we analyze the quantum solutions in the limiting cases of small and large 
scale factor. For the very small values for the scale factor the terms with coefficients g c and g\ are 
negligible in comparison with the terms with HL parameters g r and g s . In this limit the solutions to 
the equation (30) read as 



*(a) = a 



(34) 



Since the behavior of the Bessel functions for the small argument is J u (z) ~ z u + 0(z v+ ) and 
^(3-2) )i we se t c 2 — and simplify the above relation for small a as 

*(a) ~ ai+V 1 - 4 ^. (35) 

The probability of creation an universe with scale factor a is 

!a 2 , 1 - Akg s < 0, 
(36) 
a 2+ 2v /l-4feg S) 1 _ 4kgs > 0) 

which in both cases describes an expanding non-singular behavior in the early times of cosmic evolu- 
tion in agreement with the classical solutions presented in figure 2. 



8 





Figure 5: Left: Late square of the wave function (37) for k = 1. Right: The same figure for k — — 1. 

For the large scale factor, a WKB analysis like we have done for the flat case yields the classical 
solutions (27)-(29), which shows that our treatment for quantization of the model lies in a right way. 
However, for a >> 1 we can neglect the last two terms with coefficients g r and g s in (30) and write 
its solutions in terms of the Airy functions as 



*(a) = ciAi 



kg c ~ 9ka 

2/3 

. 9a 



+ c 2 Bi 



kg c ~ 9\a 

2/3 

• 9a 



(37) 



A glance at the shape of this wave function which is qualitatively plotted in figure 5 for typical values 
of the parameters, shows that it has a damping oscillatory behavior denoting a classically allowed 
region. Such a wave function is also found in [13] to describe the dynamical behavior of a universe 
dominated by the cosmological constant. 



4 Cosmological dynamics with perfect fluid 
4.1 Classical model 

Now, we assume that a perfect fluid in its Schutz's representation is coupled with gravity. In this case 
the Hamiltonian (21) describes the dynamics of the system. The advantage of using Schutz formalism 
is that, in a natural way, it can offer a time parameter in terms of dynamical variables of the perfect 
fluid. Indeed, the equations of motion for T and Pt read as 

f = {T,H} = ^, P T = {P T ,H} = 0. (38) 

A glance at the above equations shows that with choosing the gauge iV = a 3w , we shall have 

N = a 3uJ =>- T = t, (39) 

which means that variable T may play the role of time in the model. Therefore, the Friedmann 
equation H = can be written in the gauge N = a 3u as follows 

a 2 = -g c ka^ + g\a 6u,+2 + g r k 2 a 6u1 - 2 + g s ka 6 "- 4 + P a 3u) -\ (40) 

where we take Pt = Po = const, from the second equation of equation (38). Like the previous section, 
let us deal first with the solutions of this equation for the flat background, k = 0. In this case, the 
solutions can be represented by 

3v^(l-u,) 2Fl v2T^'2'2T^'-^ a j-^-^ ^ 

for 1 (for to = 1, there are no real solutions), where2-Fi (a, b, c; z) is hyper geometric function and 
to is a constant of integration. In figure 6 we have plotted the scale factor versus time for several 
values of co. As is clear from the figures the resulting cosmology is consisted of two contraction and 
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Figure 6: Qualitative behavior of the scale factor versus time based on relation (41) with numerical values Pq = 1, 
<7a = 1 and to = 1 for the parameters. The figures are plotted for oj = —1/3, —1,0, 1/3 from left to right. The dashed 
line in the figure corresponding to the case uj = — 1 shows the expectation value of the scale factor, see (59). 



expansion branches which are separated from each other by a classically forbidden region in which 
there are no physically acceptable solutions. 

For k 0, it is better to deal with the equation (40) in some special cases. To this end, we shall 
choose typical values for the equation of state parameter oj as follows. 

• uj = —1/3, for which equation (40) reads 

a 2 = -g c ka~ 2 + #a + g r k 2 a~ 4 + g s ka~ 6 + P a~ 2 . (42) 

In the region where the scale factor is small, i.e., in the early times, we can keep only the terms with 
coefficients g r and g s on the right hand side of the above equation which results the following solution 



(kg r a 2 - 2g s )^g r a 2 + kg s 

2~kg} = ±t ~ ^ (43) 

For large scale factors, on the other hand, the other terms become more important and one obtains 
the late time behavior as 

a(t) = [g 2 A (±t - t ) 2 - (P - gje)] ^ . (44) 
• uj = 1/3, for which equation (40) takes the form 

a 2 = -g c ka 2 + g A a A + g s kaT 2 + P + g r k 2 , (45) 

which has the solutions 

M 2 



a(t) = (±VgIkt - t ) , (46) 
for the early times, and 

a(t) = (±t - to)- 1 , (47) 

for late times of cosmic evolution. A quick look at these solutions shows that all of them have some 
kinds of singularities. In the next subsection we will deal with the quantization of this model to see 
how things change according to the quantum picture of the corresponding cosmology. 

4.2 Quantum model 

Now, let us to investigate how the above picture may be modified if one deals with the quantization 
of the model described by the Hamiltonian (21). The Wheeler-DeWitt equation corresponding to 
this Hamiltonian reads 



cP 1 d d ( q hp 1 o k 



*(<J,T)=0, (48) 



in which we have taken the factor ordering parameter p = 1 as before. Separation the variables in 
the above equation in the form 

*(a,T) = e iET ip(a), (49) 
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yields 



a 2 - a j- + 4 (-3c^ 4 + M^ 6 + 9rk 2 a 2 + s> s fc + Ea 3 ' 3 ^ 



da 2 da 

If k = 0, the above equation has exact solutions for some special values of u as 



V>(a) = 0. 



(50) 



V'e(a) 



Cl Ai ( -*±ffi£ 

9a 



+ C2 Bi -z±m 



273 

3a 



u = -1/3, 



C1J1/3 (fV^ + ^AO 3 ) +C 2 «/_l/3 (|V^ + 5A« 3 ) 

ClM (5y^,|.4^ a 3) +C2C/ ( 



a; = —1, 



5y/gX+2iE 5_ 4 



(51) 



w = 0, 



where M(a, b; z) and U(a, b; z) are confluent hypergeometric functions. Now the eigenfunctions of the 
Wheeler-DeWitt equation can be written as 



^ E {a 1 T) = e lET ^ E {a). 



(52) 



We may now write the general solution to the Wheeler-DeWitt equation as a superposition of its 
eigenfunctions; that is, 

tf(a,T)=/ A(E)* E (a,T)dE, (53) 
Jo 

where A(E) is a suitable weight function to construct the wave packets. Since the above relations 
seem to be too complicated to extract an analytical expression for the wave function, let us focus 
our attention on the case uj = —1 for which analytical expression for the integral (53) is found if we 
choose the function A{E) to be a quasi-Gaussian weight factor. To end this, we take c 2 = to satisfy 
the boundary condition iPe(cl = 0) = and write down the wave function as 



^(a,T) = J™ aA(£)e- i9AT e i£T J 1/3 Qv£a 3 ) d£, 
where £ = E + g\. Now, by using the equality 



z u+1 J u {Pz)dz 



-e 4 « . 



(2a)"+ 1 

we choose the weight function as A{£ ) = £ l l^e~ l£ , to obtain 

i 2 /la 6 



(54) 



(55) 



#(a,T) = Me~ i9AT 



( 7 - iT) 4 /3 



cxp 



1 7 — iT 



(56) 



where N is a numerical factor and 7 is an arbitrary positive constant. Now, having the above 
expression for the wave function of the universe, we are going to obtain the predictions for the 
behavior of the corresponding cosmological dynamics. In general, one of the most important features 
in quantum cosmology is the recovery of classical cosmology from the corresponding quantum model 
or, in other words, how can the Wheeler-DeWitt wave functions predict a classical universe. In 
this approach, one usually constructs a coherent wave packet with good asymptotic behavior in the 
minisuperspace, peaking in the vicinity of the classical trajectory. On the other hand, in an another 
approach to show the correlations between classical and quantum pattern, following the many-worlds 
interpretation of quantum mechanics [14], one may calculate the time dependence of the expectation 
value of a dynamical variable q as 



< q >-- 



< > 



(57) 
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Following this approach, we may write the expectation value for the scale factor as 

J^^*(a,T)a^(a,T)da 



< a > (T) = 



f™V*(a,T)y(a,T)da 



which upon substitution (56) yields 



<a>{T) = a 7 2 + (<7A + ^o)T 2 



(58) 



(59) 



This relation may be interpreted as the quantum counterpart of the classical solutions (41) with 
uj = — 1. However, in spite of the classical solutions, for the wave function (56), the expectation 
value (59) of a never vanishes, showing that these states are nonsingular. Indeed, in equation (59) 
T varies from — oo to +oo, and any To is just a specific moment without any particular physical 
meaning like big-bang singularity. Now, let us take a look at the uj = — 1 case of the figure 6 in which 
the expectation value (59) is plotted with the dashed line. As is clear from this figure, for a perfect 
fluid with co = — 1, the corresponding classical cosmology admits two separate solutions, which are 
disconnected from each other by a classically forbidden region. One of these solutions represents 
a contracting universe ending in a singularity while another describes an expanding universe which 
begins its evolution with a big-bang singularity. On the other hand, the evolution of the scale factor 
based on the quantum mechanical considerations shows a bouncing behavior in which the universe 
bounces from a contraction epoch to a reexpansion era. Indeed, the classically forbidden region is 
where the quantum bounce has occurred. We see that in the late time of cosmic evolution in which 
the quantum effects are negligible, these two behaviors coincide with each other. This means that the 
quantum structure which we have constructed has a good correlation with its classical counterpart. 

For a background geometry with k 7^ 0, to analyze the quantum behavior of the model we may 
neglect the terms with coefficients g c , g\ and g r in (50) in the early times, i.e., in the region where 
the quantum effects have their dominate role. In a such a situation this equation takes the form 



da 2 



ld_ 

a da 



+ 4 



9sk 



+ Ea 



l-3u> 



Ha) = 0, 



(60) 



where its solution for some special cases are as follows 



(61) 



in which we have again applied the boundary condition ip(a = 0) = on the eigenfunctions. Following 
the same steps which led us to the wave function (56), we obtain the wave function as 



„l+\A-4fc 9s 



( 7 -iT) 



4( 7 -iT) 



-1/3, 



ji+V 1 - 4 *^ 



• exp 



*(a,T) = < 

k ( 7 _j T )l + V 1 - 4fc 9s 

from which the expectation values are obtained as 

r ( 7 2 +t 2 ) i/4 , 



(62) 



7-iT 



w = 1/3, 



< a > (T) ~ < 



•1/3, 

^ ( 7 2 + T 2)l/2 5 w= l/ 3 . 



(63) 



The discussions on the comparison between quantum cosmological solutions and their corresponding 
form from the classical formalism, i.e., equations (42)-(47) are the same as previous model, namely 
the flat model. Similar discussion as above would be applicable to this case as well. 
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5 Summary 



In this paper we have applied the recently proposed Hofava theory of gravity to a FRW cosmological 
model. After a very brief review of HL theory of gravity, we have considered a FRW cosmological set- 
ting in the framework of the projectable HL gravity without detailed balance condition and presented 
its Hamiltonian in terms of the minisuperspace variables, both for the vacuum and perfect fluid cases. 
For the flat model without the matter contribution, we showed that the classical field equations admit 
contracting and expanding de Sitter-like solutions in which the cosmological constant is modified by 
the HL parameter A. For the non-flat background in this case, though the corresponding Friedmann 
equation did not have exact solutions, we analyzed the behavior of its solutions in the limiting cases 
of the early and late times of cosmic evolution and obtained analytical expressions for the scale factor 
in these regions. We saw that these solutions are consisted of two separate branches each of which 
exhibit some kinds of classical singularities. We then have repeated the calculations when the model 
is augmented with a perfect fluid as the matter field. Again, we showed that the classical solutions 
have either contracting or expanding branches which are disconnected from each other by some clas- 
sically forbidden regions. Another part of the paper is devoted to the quantization of the model 
described above. For an empty universe, we have shown that by applying the WKB approximation 
on the Wheeler-DeWitt equation, one can recover the late time behavior of the classical solutions. 
For the early universe, we obtained oscillatory quantum states free of classical singularities by which 
two branches of classical solutions may communicate with each other. In the presence of matter, we 
focused our attention on the approximate analytical solutions to the Wheeler-DeWitt equation in 
the domain of small scale factor, i.e. in the region which the quantum cosmology is expected to be 
dominant. Using Schutz's representation for the perfect fluid, under a particular gauge choice, we led 
to the identification of a time parameter which allowed us to study the time evolution of the resulting 
wave function. Investigation of the expectation value of the scale factor shows a bouncing behavior 
near the classical singularity. In addition to singularity avoidance, the appearance of bounce in the 
quantum model is also interesting in its nature due to prediction of a minimal size for the correspond- 
ing universe. It is well-known that the idea of existence of a minimal length in nature is supported 
by almost all candidates of quantum gravity. 
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